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&: ABSTRACT 



We present a new study of polarized elastic niuon-electron scattering. The Born cross-section is 

• 5^ \ calculated for arbitrary polarization of muon and electron. The complete photonic 0{a) radiative 

/\ ' corrections are determined for the case of longitudinally polarized muons and electrons. All calcula- 

C^ ■ tions are done by two methods: semianalytic, which allows an implementation of the experimental 

cuts used for the analysis of //e scattering data from the beam polarimeter of the SMC experiment 

at CERN and completely analytic, which is used for cross checks. The FORTRAN code fiela realizes 

formulae of both approaches. We prove that certain experimental cuts lead to negligible radiative 

corrections in the muon beam polarization experiment. 
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1 Introduction 

Polarized elastic /le scattering is being measured by the SMC collaboration at CERN as a 
monitor of muon beam polarization [||]. Since the measurement pretends to be very precise, 
the photonic corrections have to be taken into account. 

The differential cross-section for this process in lowest order may be cast into the simple 
form 101 



, BORN _ g 

da Zrca 



^^ (1 - yP.P,) + 1(1- PeP.) 



:i.ii 



dy rrieEf, ,_ 

where the following notation is used: 
m^, iTLf. - muon and electron masses, 
P^, Pe - longitudinal polarizations of muon beam and electron target, 

K 

y = y^ = I — -7- - the measured energy loss of the muon, 

(Tfi \ 
1 H — - its kinematical maximum, 

E^, E , E^ - muon (initial, final), electron final energies in the laboratory frame. 
The polarization dependence of da is used to calculate the measured electron spin-flip 
asymmetry AJ^J^ 

dajU) c?a(TT) 

/lexp _ dy dy ,-. „% 

^'^^ da{]i) ^ dam • ^ ' 

dy dy 

The asymmetry is measured as a function of the variable y^. 

Previous calculations @]-[^, presented results in terms of the variable ye = E^/E^, and, 
only Ref. took into account the polarizations. 

If elastic fie scattering is treated in the Born approximation then 

yf, = ye = y, (1-3) 

and eq. (|1 . 1| ) may be written in terms of either y^ or y^. 

The situation changes drastically if one calculates QED corrections. Due to the emission 
of non-observed photons the identity (|1.3| ) does not hold anymore, and one has to specify the 
variable to be used for the calculation of radiative corrections. Their numerical values may be 
very different in y^ and ye- 

Since the measurement and the analysis were performed in terms of y^ |l| , the calculation 
of QED corrections must be done, of course, in terms of the same variable. This is why a new 
calculation was neccessary. 

Our new calculation is the theoretical basis for the Fortran program fiela, [0. It is a 
complete, order 0{a^), calculation. It takes into account longitudinal polarizations of both 
/i and e, finite muon mass effects (the electron mass is neglected wherever possible). In the 
semianalytic approach it is possible to apply all experimental cuts which were used in the 
analysis of the experimental data: 

- a recoil electron energy cut, E^> E {E =35 GeV); 

- an energy balance cut, E — E^ — E^ > E {E =40 GeV); 



angular cuts on both /i and e, ^^ and 6^ in the laboratory system: 

BORN 

< H ■ i H ■ = I mrarl I In thp a hnvp H 



o: 



e: 



< 



/imcas n 



while 9,, and 9 



< ^TT,i„ f^min = 1 mrad). In the above, 9'^'^^ and ^^'^^^ are the measured angles, 
are angular values calculated using BORN kinematics. 



min V iiiin 
BORN 



Muonic RC 




,+ 



+ 



+ 



Electronic RC 




Vacuum Polarization 




> 



Figure 1: Feynman diagrams for the elastic /xe scattering in order 0{a^) 
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It order 0{a^), the 14 Feynman graphs, shown in Fig. |T], contribute to the cross-section. 
The latter may be subdivided into 12=2x6 separately gauge invariant contributions: 

7 QED 2 6 J I 

"i/M 1=1 k=l "i/M 

where the indices k and / have the following meaning 

1=1 — unpolarized contribution, / = unpol; 

2 — polarized contribution (the terms proportional to PeP^ in ( |1.1| )) , / = pol. 

k = 1 — Born cross-section, k = b; 

2 — Radiative corrections (RC) for the muonic current: vertex + bremsstrahlung, k = nn; 

3 — contribution of the anomalous magnetic moment of the muon, k = amm; 

4 — RC for the electronic current: vertex + bremsstrahlung, k = ee; 

5 — fie interference: two-photon exchange -|- 

muon-electron bremsstrahlung interference, k = fie; 

6 — Vacuum polarization correction, running a, k = vp. 

The resulting QED corrected cross-section is given by the sum 

J QED / J unpol 7 pol \ 

dj^^Y. P^ + P.P.^]- (1.5) 

The cross-sections with k = fifi, ee, fie have similar generic structure 

, , BORN , BREM 

dak a vRda da^ ,. „. 

-^ = -^k -1 + -; . (1-6) 

dyf, vr dy^ dy^ 

with a factorized part 6^^ originating from infrared divergent virtual (V) and real soft photon 
(R) contributions^. 

The main conclusion of this study is illustrated in Figs. ^ and |^, which present radiative 
corrections to the asymmetry as a function of the variable y^ for two cases: without any cuts 
and with experimental cuts described above. 

The asymmetry A" , and the radiative correction to it, 6y are defined as follows: 



da'^iU) da'^i]]) 



^^'^ " ^%I) i%Ty' "" " ^°^^' ^™' ^'^ ~ A^' 



QED 



dya dy^ 



fie 



As is seen from the figures, the corrections without cuts are very large and reach up to 
-20%. When the four above mentioned cuts are taken into account they reduce S to values 
below 1%. Actually, for a wide range of |/^ they are even well below 1%. 

The main conclusion of our investigation is that one may safely neglect radiative corrections 
in the determination of the muon beam polarization with the SMC set-up. 



^ In the following we will always present the formulae in the form (1^), i.e. summed over I 



u 

-2 


- 0%] 










ju-ela 


-4 


- 












-6 


-\ 












-8 
10 


\ \ 












14 
16 


\ 


^ 


\ 


^ 




J 


18 
on 


\ 


1 1 1 1 1 


1 1 1 1 1 


1 1 1 1 1 


1 1 1 1 1 




-% 


1 0.2 0.3 


0.4 


0.5 


0.6 


0.7 


0.8 0.9 1 



Figure 2: QED corrections to the polarization asymmetry without experimental cuts. 
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Figure 3: QED corrections to the polarization asymmetry with experimental cuts: 
E = 35 GeV, E =40, ^e,min = ^^. min = 1 mrad. 



2 The Born cross-section 

2.1 Kinematics and phase space 

We consider the elastic scattering process 

fi{h) + e{pi)^fi{k2) + e{p2), (2.8) 

in a fixed target experiment, i.e. with the initial state electron at rest, ki = (0, ime). 

Since the typical incident muon energy, E^^, in present-day fixed target experiments is 
0(10^ — 10^ GeV), the maximal c.m.s. energy, 

s = -{ki+pi)^ =ml + ml + 2meE^, (2.9) 

is very small, y/s < 1 GeV. Therefore, we may completely neglect Z-boson exchange. 

In fact, for the energy used by the SMC collaboration, E^ = 190 GeV, the invariant s is 
only 20 times bigger than the muon mass squared. Therefore, we can not neglect effects of 
the finite muon mass. Of course, the electron mass may be completely neglected. 

While calculating the Born cross-section, we will perform all derivations exactly even in 
me, since the resulting expressions are very compact even if rUe is kept, but at the end of 
calculations we will neglect the electron mas^. The Born process is characterized by one 
kinematical variable, besides s. We will use the dimensionless variable y: 

2.10) 



2.11) 

2.12) 
2.13) 

2.14) 

2.15) 
2.16) 

2.17) 



Pi{ki - fca) 




- 1 


K 


Pih 


k^ " 




E,. 



y 

We will introduce also the transferred momentum squared 

It is easy to derive the identity: 

Q' = Sy, 
where 



,2 _2 
= S — If 

For the Born cross-section, we have 



O 2 2 

b = s — m^ — m^. 



with 



As 



X^ = S^ - 4mlml, 






In terms of y the differential phase space reads 



dT2 = —r^^=d(^ = j=dydLp. 



We refer to this approximation as to the "Uhra-Relativistic Approximation (URA) in me" ■ 



2.2 Spin degrees of freedom 

Since we are going to deal with the scattering of polarized particles, there will be additional 
essential variables, besides s and y, which are supposed to describe the spin degrees of freedom 
of the problem. Their description uses the language of spin density matrix ( for details we 
refer e.g. to Appendix C of 0). 

For non-polarized particles, we use projection operators in trace calculations, i.e. summing 
and averaging over spin indices looks as 

Y.u%p)u%p) = Uip), (2.18) 

with 

A{p) = -ip + m. (2.19) 

For polarized particles, we use the spin density matrix instead 

Y.u%p)u%p) = 1(1 + ^750A(P), (2.20) 

s ^ 

where ^ is the polarization four- vector 0. 
In the particle rest frame, p = 0, it is: 

e=(Pn,0), (2.22) 

where n is a unit vector in the direction of spin quantization, and P is the polarization^ 
defining the degree of spin orientation along the direction n. For instance, P = 1 means that 
the probability of a particle to have its spin projection along the direction n is equal to 1 (right 
handed longitudinal polarization, if vector n is chosen along particle momentum p\ From 
( p.22| ) in the particle rest frame, we have 

ip = 0, 

^2 = p2_ (^2.23) 

Due to Lorenz invariance, the properties (|2.23|) are fulfilled in any Lorenz frame. 

The initial electron with the four-momentum pi is at rest in the laboratory frame. Using 
then the direction of incoming muon as the direction of spin quantization, i.e. n = ki/\ki\, we 
get the four- vector of the electron polarization from ( |2.22| ) 




ee = Pe -^,0 . (2.24) 



The four-vector ^^ may be obtained from the expression similar to ( p.24| ) in the muon rest 
frame 

Cm = ^.1^,0 I (2.25) 



^ A naive use of longitudinal polarizations from the early beginning of calculations, i.e. use of the spin 
density matrices in the form 

J2u%p)u%p) = ^{l + Xj,)A{p), (2.21) 

S 

does not properly reproduce the finite muon mass terms in the PgP^ part of the cross-section. 



by Lorenz boost to the electron rest frame along the beam axis: 

A;? / ki \ki\ 



^p 



P„ 



m„ 



Ifcll' fc? 



(2.26) 



We will consider the Born cross-section with arbitrary orientations of electron and muon 
spins. We choose the laboratory frame with 2;-axis oriented along the incoming muon 3- 
momentum ki and with the plane {x,z) coinciding with the reaction plane. Another plane is 
spanned by the vectors ki and the projection {^e)xy of vector ^g to a plane perpendicular to 
2;-axis. In this frame, the relevant 4-vectors are written as follows: 

k, = (O, 0, |fci|,fc?), 

h = (|^2|sin6'^, 0, |A;2| cos6'^, /c^j , 

pi = (0, 0, 0, me) , 

P2 = (IP2I sin 616, 0, |p2 1 cos 6*6,^2) • (2-27) 

For the spin vector C,e arbitrarily oriented in 3D-space, we have in the choosen laboratory 
frame instead of (|2.24| ) the following generalization 

^e = -Pe (sin-i^eCosyje, siu "i^e siu (/?£, cos-f^g, 0). (2.28) 



We may identify the angle cpe in ( |2.27| ) with ip of the phase space parametrization in ( |2.17| ) 
and therefore if becomes an essential degree of freedom in presence of a transverse polarization. 
For arbitrarily oriented ^^ , we have instead of ( ^.25| ) in the corresponding rest frame 



^,1 = Pfi (sin ^^ cos ipf,, sin t9^ sin yj^, cos ^f,, 0) 
Now we boost ^„ from the muon rest frame to the laboratory frame 



e. 



P^ I sin -(9^ cos (^^, siw&f^smip^, ^-cos^ 



\ki\ 



m,, 



cos -(9, 



(2.29) 



^2.30) 



Using the explicit representations ( |2.27| ), ( p.28| ) and ( |2.3CI| ), we can easily derive all scalar 
products involving polarization vectors. 

The doubly differential (in y and cp fl) cross-section exact in both masses reads: 



da 



dyd^ 



2a'^S [1 si 



1 

y 



1 

Y 



K \f ys ' 2 ' ^^^'^ 

] ]=^ I 1 H ^ smt^e cos -(76 sm-fy^cos*/?^ 



1 H — =■ ] sin 9^ cos f?^ sin t9e cos 99 ^ 



- I cos df. cos d^ 



yy^s 

me\k2\ 



y\^ 



-2 — r--^ sm ■ii'e sm i^^ 



N|P2 

s 



sin 9e sin 9^ cos ipe cos ip^j, 



+2/ cos 5(^^_^^) 



(2.31) 



For ip two choices are possible: 1) ip — ipe, then i/j^ = (pe^ S,ip^-ip ,) or 2) p ^ p^, then p^ = <P/^+ 5 



CPc-ffi)- 



The expressions for sin O^, and sin 9^ exact in irtf, are 



sin6'e 
sin 0,,, 



2mey/Sy 



(2.32) 
(2.33) 



where 



and 



A/ = 


- S\l-yr~Amlml 


A° = 


= SV + ^rnlSy, 


y = 


- 'V'y) 


Y -- 


sS \ S 



(2.34) 
(2.35) 



(2.36) 



is the kinematical maximum of y- variation. 

The substitution of these variables into ( 2.31|) exhibits an interesting property of the general 
Born cross-section which becomes 



da 



dydip 



2a^S\l__^ 1 



A. 



y^ 



ys 



PeP^. 



h — - - COS ^e COS ^, 

. y y ^ 



ni^y/Sy 



y\JK 

rriey/Sy 



y\^s 



,!+¥ 
\.'^ 



COS "Oe sin 'd^ cos ip^ 
cos 'd^ sin -de cos ipe 



—2 ^ sm Vp sm v. 



yS 



1- 



Y 



cos ipe COS (p^ + COS 5(<^^_^,^) 



(2.37) 



From the last presentation it is clearly seen that while terms related to the transverse electron 
polarization are small since they are suppressed by the electron mass (third and fourth lines), 
the term induced by the transverse muon beam polarization (second line) is not small, since 
it appears to be proportional to the muon mass. 



By trivial algebra one may show that the expression (|2.31|) reduces to a very short form in 
two particular cases. In case of tranverse electron and longitudinal muon polarizations in the 
URA in the electron mass we obtain 



da 



dydip 



2«2 



111 ^ 1-y fl I 

y^ yY 2 y \2 Y 



(2.38) 



The corresponding expression for the case of longitudinal polarization of both particles, but 
exact in both masses reads: 



da 



BORN 



dy 



A. 



4-^ + ^ + PeP, 



yS 



e-^ fi 



1 
y 



1 

Y 



V 



(2.39) 
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Having in mind that r^ = a /me and S = 2mf.E^, we immediately identify (|2.39| ) with the 
corresponding expressions from Section 2.2 of ref. if one neglects terms of 0{ml) here. 
In the URA in me, equation ( p.39|) may be rewritten in a form, which is explicitly positive 



definite: 

BORN 



da 4:Tra^ 



dy S 



(y-y) ,. ..r^r.^ , 1 



l-yPM + -[l-PeP. 



y2Y ' . e M, , 2 



(2.40) 



3 Complete 0{a) Radiative Corrections 

3.1 Kinematics of fie -^ fie^ 

The reaction 

fi{ki) + e{pi) -^ fi{k2) + e{p2) (3.1) 

is accompanied by the bremsstrahlung of non-observed photon(s) 

/i(fci) + e(pi) -^ /i(/c2) + e(p2) + {n)-i{p). (3.2) 

First of all we will study the kinematics of one-photon bremsstrahlung. We want y^ to be the 
last integration variable out of a set of four variables (besides S). There is some freedom in 
doing this. 



(3.3) 



We will use the definitions: 




Ql = {k, - k2)\ 


Pi{ki-k2) 
Piki 


and 




Ql = {P2-Pl)\ 


^, PliP2-Pl) 

Pih 


The other invariants are: 




zi = -2pki, 


Z2 = -2pk2, 


Vi = -2ppi, 


V2 = -2pp2. 



(3.4) 



(3.5) 

Of course, only four of the invariants are independent. 

Using 4-momentum conservation, it is easy to derive the following relations among them: 



zi + V^ 


= Z2 + V2, 


Vi 


= Sy.-Ql 


V2 


= Sy,-Ql, 


Ql 


= Sy,. 



(3.6) 
We will use the following set of independent variables 

S, Q% y„ Ql ^2(1). (3.7) 



The last of eqs. ( |3.6| ) deserves a comment. Our choice of the 4-momentum pi in defini- 
tions (|373|) and ( |3.4|) introduces the asjTiimetry between y^ and ye- This is reason why Q^ and 
y^ may be chosen as independent variables, while there exists a relation between Ql and ye- 



ll 



For the moduli of particle momenta and the energies, in the electron rest system pi = 0, 
we obtain 

\p\ = ^, / = ^^y^~y'\ 



I "'I I ~ T, 5 ""l 



2mp 2mp 



''e 



I r I V '^^ 7 "^i 



2mp 2mp 






(3i 



2mp ^ 2m, 



e 



■)2 



IP2I = |<5e| = 1^ , P2 = "^e + 



^5i Q,. 2 r - S{Ql + 2ml) S-QI-Z2 



ki ■ ko = n — m,,, ki ■ P2 — o 

with 



4m^ 2 ''' " 4m: 



^i = 5(l-y^). (3.9) 

The corresponding equations for the Born kinematics (|p| =0) can be easily derived from 
(|]8D setting 

Vi = V2 = z^ = Z2 = Q. (3.10) 

In this limit 

Ql = Ql = Q'' and ye-y>. = y. (3.ii) 

Many relations and useful notations can now be taken from [§. As usual, we introduce 
the relevant kinematical A-functions: 

A, ^ X[-{p^+p)^-pl-p'] = S\y,-yp)\ 

A, ^ Xl-^p, + k^)\-pl-kl] = S'-4mlml 

A, = X[-ip^ + k2)^-pl-k|] = S!-Amlml (3.12) 

A^ ^ X[-{p, + Q^f,-pl-Ql] = S'y.' + AmlQl, 

Xe = \[-{pi+P2y,-pl-pl] = SW + AmlQl 



where 



\{x,y,z) = x^ + y^ + z^ -2xy-2xz-2yz. (3.13) 



3.2 Kinematic boundaries 

The boundary conditions may be taken from 0. The first one, (B.3) of 0, remains unchanged 

mlQt + S\Ql + mlS\' - XsQl = 0, (3.14) 
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while the second condition, (B.4) of 0, changes, since, contrary to 0, ^/e is not an independent 
variable here. This is due to the fact that we are dealing now with elastic scattering rather 
than with the deep inelastic scattering in 0. Using the last of eqs. (|3.6|) , eq. (B.4) of takes 
the form: 

S%'Ql + QtQl - mliQl - Qlf - Sy,Ql{Ql + Ql) = 0. (3.15) 

The physical region £^ = {Q'^^.y^) is given by two inequahties (see 0, subsection B.2.1), 
which are derived from ( |3.14| ) 



A. 



Q<Ql< o , J , ^2 ^ Ql- (3-16) 



'^' ~ S + ml + ml 



where 

Q 

S 



vr^QD = ^ <y,<yT''iQl)^ (3-i7) 



vrm) - ^ 



^yA^Am - Ql 



(3.18) 



and 

An. = Ql{Ql + ^ml). (3.19) 

The solution of eq.( p.l5D is 



(g2)max,min ^ ^^^ '^ ^^^J^^ U^ \' '^ ' ■ (3.20) 



Sy,{Sy, - Ql) + 2mlQl ± jS y, - QD^X 
2{Sy,-Ql + m: 



3.3 Another set of independent variables 

Besides of ( p.7| ) we use 

S, y^, V2, V^i, ^2(1), (3.21) 

To write down the limits in these invariants, we reorder first the physical region £^^ = 
(Q^)1/m) ~* iVfi^Qfi)- Trivial manipulations with ( p.l6| )- p.l8| ) lead to 

0<y,< l/;^^ (3.22) 

iQir^' <Ql< min{(g^)--, Sy,}. (3.23) 

Here 



A., — S'^yn ± a/ A„A; 
''^^ 2mi 



(g2)max,mm _ ^ ^^2 (3-24) 

Solving the equation 

iQT^'' = Sy„ (3.25) 

we find a maximal value y™^^ 

C = 7777^ ^^ 9T, (3-26) 
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where the two upper hmit branches of ( 3.23|) meet each other. 

From ( |3.23| ) and the definitions of Vi, we easily derive the hmits of V2 as function of y^ 



0< Vo < 



1ml 



(3.27) 



The second solution 



V^ 



2ml 



(3.28) 



is unphysical. It is negative in the physical region (|3.22| ) of y^. Finally, from ( p.20| ) we derive 
the limits of Vi as functions of y^ and V2 



^r'° <Vi<V, 



1 5 



where 



VI 



max,min -rr 

— V2 



Sy^ + 2ml±J\, 



2{V2 + ml) 
Examining (|3.27| ), one may see that the invariant V2 is positive only in the interval 



^<y,<yl 



(3.29) 



(3.30) 



(3.31) 



To complete the study of kinematics of the reaction ( |3.2| ), we have to give the limits of 
variation of the variable -21(2)- We may take all the relevant formulae from |^ and simply list 
them for completenes: 



^1(2) li/M' V/X5 iJe) 



5i(2) ± /D; 



Cl(2) 



^1(2) 5i(2) T ^/Dl 

and the Gram determinant 

R, = -Aizf + 2BiZi-Ci = -A2zl + 2B2Z2-C2, 
A2 = Xf, = Ai, 

B2 = {2mlQl{Ql-Ql) + S\l-y,){y,Ql-y,Ql) + S\l)Ql{y,-y,)} 
^ - B,{il)^-il-y,)}, 



(3.32) 



(3.33) 
(3.34) 

(3.35) 



C2 = {Sil-y,)Ql-SQl[il)-y, 
= C,{{l)^-{l-y,)}, 



+ 4m? 



D. 



Bl{2) -^l(2)C'l(2)- 



S\y, - ye){yM - yeQD - ml{Ql - Qlf) 

(3.36) 
(3.37) 



Here we now understand that Ql = Sye- This makes no simplifications, therefore we did not 
substitute it in order not to destroy a nice symmetry of these equations. 

The inequalities ( |3.31|) , (|3.27|) , ( |3.30|) and (|3.32|) are the kinematical limits in the sequence: 

1//., V2, Vi, Z2. 
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3.4 Phase spaces 



In order to have a convincing proof that all the boundaries are correctly derived, we always 
write a supporting FORTRAN program to numerically check that the phase space volume in 
any sequence of variables is exactly the same. We calculated as many phase space integrals 
analytically as is possible. 

For example, we realized the cross check of our sequence ( p.21|) . Three integrations may 
be easily performed analytically, yielding: 



r 









max _ ^2 1^ _^ e 



m" 



(3.38) 



The numerical check returned for this integral exactly the same value as for the basic sequence, 

(c.i)ofi. 

We will always integrate over Z2 first. Then, only three variables Vi,V2,y^ remain. The 
physical regions £ = {y^, V2) and X = (V2, Vi), as derived from (|3.27| ) and (|3.30|) , are shown 
in Figs. 4 and 5. 





Figure 4: {y2, y ^i) -plot 



Figure 5: {Vi,V2)-plot 



3.5 Bremsstrahlung cross section 

For the normalized bremsstrahlung cross-section, we have 



da 



2JX 



M 



dT, 



with the 2^3 phase space 



dTs = (27r) 



d^k2 d^p2 d^p 

{27i)^2k^2 WV2 (27r)32p' 



:S{ki+pi -k2-p2 -p)- 



(3.39) 



(3.40) 
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In terms of the variables ( 3.21 ) the differential phase space looks as follows 



dV; = ^dy,dV2 dVi^. (3.41) 



We developped two branches to present the final result: 

• calculation within the Numerical Approach, with the possibility to impose arbitrary 
experimental cuts; 

• calculation within an Analytic Approach with limited possibility to apply cuts. 

We note that the integrals over zi were calculated analytically with arbitrary limits ;g™™'™'^^ 
(see, Appendix (^), while integrals over V2(i) were computed purely numerically. This is the 
essence of our numerical approach, the muenum branch of fiela. 

We also performed the complete 0{a) RC calculations within an Analytic Approach, the 
mueana branch of fiela. In this description, we will present only a sketch of the calculations. 

Within the analytic approach we have integrated the bremsstrahlung contribution over the 
variables Zi and Vi over the full photonic phase space without imposing experimental cuts. The 
aim of this is twofold: first, we performed an independent calculation of the bremsstrahlung 
cross section and its integration which was used to cross check the formulae of the numerical 
approach and their coding in the FORTRAN program fiela. Second, these analytic formulae are 
rather elegant, although lengthy, in the Single Ultra Relativistic Approximation (SIURA) in 
the electron mass. Due to analytic integration, they run at the computer incredibly fast, several 
orders of magnitude faster than the numerical integration within the numerical approach. We 
note that this relative simplicity may be reached only if we integrate over the full phase space 
of photons. From formulae in the SIURA it is easy to derive the formulae in the DOubly 
Ultra Relativistic Approximation, DOURA. It is important to stress, that in the latter case 
satisfactory precision of the calculations is not ensured. 

The derivation of analytic results is quite straightforward. We consequently integrate in the 
sequence dV2dVidzi. In reality, it is convenient to calculate and substitute the twofold integrals: 
dVidzi. These are presented in the Table in Appendix (|D.1|) . After substituting the dVidzi 



integrals, we arrive at another Table of integrals over V2 which are listed in Appendix ( p.2|) . 



They are much more complicated than the first, two fold integrals. The latter integrals were 
calculated in limits [0, V'g™'^^] with an arbitrary upper limit. This allows us to impose V2- 
dependent experimental cuts even within the analytic approach. An example is the muon 
angular cut. 

A comment on DOURA is in due here. Since here we neglect m^ , a problem might arise 

from the fact that Q^ is proportional to the muon mass. Apparently, m? can't be neglected 

in Qu™™) since it enters in the denominator of the photon propagator. However, we assume 

that a cut is imposed on V2, which prevents Q^ from reaching its kinematical limit. In 
other words we assume that 

Ql > iQir' = (Qir » ml (3.42) 

3.5.1 Final expression for the normalized bremsstrahlung cross section 

Defining the z-integrated contributions. 



, ,BREM 

M 



dzi 



^BREM 



(3.43) 
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and neglecting rrie wherever possible, we finally obtain 

BREM ^ „T>max _T>max 



da 



r\ „T/max 

a r2 



dVfi 



S Jo 



dVo 



V," 



ym 



dVi 



M 



(3.44) 



■^N^, ^™"^'™^^ QjQ functions of experimental cuts. In this way, the 



In (§^, the hmits V^{ 
cuts are implemented within our numerical approach. 

There are three gauge-invariant contributions to the bremsstrahlung differential cross sec- 
tion: muonic, electronic and fic interference, each of them is represented as the sum of an 
InfraRed divergent contribution 'IR' and a finite (Regular) part 'R': 



(3.45) 
(3.46) 
(3.47) 





, ,BREM 

M 


2 
lifi 


z 


= 


^l 


( 




+ 5 


R\ 
iifi) 5 








M 


2 
ee 


z 


= Ql{BFlf + Sf:), 




, ,BREM 

M 


2 
fie 




= QfiQe{BFl^ + Sl), 


where ,B is the Born factor 






y2Y 


(1 yP.P,) + 2{l PeP,)-^^l,^, 


M 


The Fj are infrared fac 


;ors: 






{Sy, + 2m^^ 


1 

Z1Z2 


Z 


" 1" 

A. 


z 


"1 

A 


5 

z 


pIR 


( Sy, ml 
[V^V2 V^ 


ml\ 

vi) 


1 


1 

z 


r.IR _ 

^ fie 


S 


' 1 




\ 


S 

^2 


' 1' 

Z2 


Z 


+ 


Si 


"1" 

_^2_ 




"1" 





(3.48) 

(3.49) 
(3.50) 
(3.51) 



In (|3.49| )- (|3.51| ) and in the cumbersome functions S^ given in Appendix (^), the z- 



integration ( |3.43| ) with arbitrary limits is assumed to be done 
Tables of 2;-integrals with cuts and equations for limits V^^ 
mental cuts are presented below in Appendix B: (§]^), (§]^), ( |B.1U71 )-( |BTI?BD 



Tables of 2;-integrals with cuts and equations for limits V^^f^ ^mm,max ^^ terms of experi- 



3.5.2 Treatment of the infrared divergent part 

The three terms with F^^ in (|3.45| )- (|3.47| ) cannot be simply integrated in ( |3.44|) because of 
the infrared divergency at V2 = 0. It is treated by dimensional regularization. 

Substituting all terms with F^^ into (|3.41| ), we define the IR Part of the bremsstrahlung 
cross-section. 



da 



IR 






^ {qXI + ^^^^K^ + Qe^if ) dT, 



a 



5" 



^ (QlKl + QMQe^i^ + QlFlf) \e{s - p') 



= da^^''°^^ + da^^'^''"^. 



e{p^ - e] 



dV. 



(3.52) 
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We will treat da^^ in the R-fraine which is defined in Appendix by the condition 

P2+p = 0. (3.53) 

In this frame we find 

V2 = -2p.p2 = -{p2 + pf -ml = (pI + p^f - ml (3.54) 

For sufficiently small e we have at the point of separation of soft and hard photons, p^ = e, 

^2 = ^2 = 2mee, (3.55) 

which can be choosen to be much smaller then any typical invariant of the process. 

From a study of the R-system it may be understood that there is an unique correspondence 
between angular integrations in the R-system and invariant integrations, i.e. we may use 
invariant limits of z in (|3.43| ) and of Vi in ( p.44| ) instead of angular limits in the R-frame in 
order to compute the hard part of da^^ 



da 



IR.hard 



7^3^,3 






B (QlFll + Q.QeC + ^l^ef) dTsOip' - e). 



(3.56) 



This is very convenient, since hard photons are in general affected by experimental cuts and 
the cuts are implemented in our approach by means of limits of a numerical integration over 
zi and Vi. Finally, we get from (|3.52|) 



^^IR,hard 

dy,. 



b Jv2 

, BORN 

da a 



dVo 



V7 



ymm 



dVi q2 






+ QuQ 



^He 



7IR 
fie 



+ Q 



:.IR 



■IR,hard 



dy^ 



TT 



with 



s:IR,hard 



dV. 



V2 



yrr 



dVi q2 



7IR 
ij.fl 



+ Q«Q 



/j"ie 



7IR 
fie 



+ Q 



:iIR 



(3.57) 



(3.5^ 



Again, in ( |3.57| )-( p.5^ ) the ^-integration ( |3.43D with experimental cuts is assumed to have been 
done. 



3.5.3 A short form of the completely differential bremsstrahlung cross-section 



The IR finite 2;-integrated contributions ( [A.71[ - |/Ur^ ) are rather cumbersome indeed. For the 
unpolarized squared matrix elements of (|3.45| - |3.47|) with not yet separated infrared parts (p. 49 



3.51 ) a compact representation is known in the literature |ni|. We present it here, in notations 
of our paper, even in a more compact form: 

- unpol 



BREM 
^form 


2 
MM 


, ^BREM 
^form 


2 
ee 


, ,BREM 
^form 


2 



Q 



^2rT- + 7^H^MM + ^eA 



QiziZ2 Qt 



2\e 
lifj. 



unpol 



unpol 



T 2 

u _ fm^A^ + m^ A'' 



QiViV2 Q, 



QmQ. 



/i"ie 



r 



s 



S' 



u 



U' 



QlQl V^l^l ^2^2 V1^2 ^1^2. 



2ml 



iQr'\ 



QiQ 



,(3.59) 



where 



and 



r 






A'' 



A'' 

ee 



*/ie 



52 + 5/2 ^ f;2 ^ f;/2 
5'2 ^ f;2 52 ^ f;/2 



^2 



^2 
^2 



-2 



-2Q^ ( - - - 

.Z\ Z2 



SU + S'U' 

Z\Z2 



Zl 



Z2J 



S 



U' S' 

+ — 



u 



Zl 



+ 4 



^IQI 



Z2 Zl Z2 

S^ + U^ 



V? 



\/2 

Vi V2 

{S + Uf + {S' + U'f 



+ 2miQ 



2 



y? 



V1V2 



+2- 






s 



S' 



Zl 



\ViZi 
U' -S 

Z2 



+ 



u 



+ 



V^2l 



U' 



V2Z2 V1Z2 Z1V2 



(3.60) 



S' -- 


= S-zi-Vi, 


u = 


= -S + Ql + Vi 


U' = 


= -S + QI + Z2. 



(3.61) 



The representation ( |3.6CI|) - (|3.61|) is really very elegant, and it may be used in a Monte Carlo 
code. It can't be used, however, within the numerical approach of this paper for the following 
reasons: 

1. It does not take into account polarizations (main reason). 

2. It includes the IRD part which has to be separated out. 

3. Our bremsstrahlung cross-section, as we emphazed already, is integrated once over the 
invariant 2:1,2 within arbitrary limits (by making use of a table of indefinite integrals). 
To do so, we need the canonical representation 



S!" = J2 fiizj)KiiiS, y,,, V2, Vi). i = ee, e/x, /i/i. 



1,2, 



/ 



24, 



(3.62) 



see Appendix (0). We gain at least one order of magnitude in CPU-time, since instead 
of a three-fold numerical integration over V2, Vi and 2:1,2, we need only two- fold: (V2, Vi). 
This gain of CPU-time makes our code very fast and user friendly. In order to arrive 
at this level, however, one has to substitute S", f/, [/', and the compact expressions, 
actually, become much more cumbersome (see Appendix A). 



4. Finally, the canonical representation ( |3.62|) seems to be the only road towards an ana- 
lytically integrated differential cross-section, see the discussion at the beginning of this 
section. 
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3.6 The muon anomalous magnetic moment contribution 



The contribution of the anomalous magnetic moment of the muon to the cross-section can be 
expressed as (here we follow decomposition of ref. |TT|] ): 



where 



and 



da' 



dy 



pamm 

s ^ 



T-'amm 



' Sy\ 



^ - 7) (2 - yPeP.) , 



1, /3 + 1 



^ 



\ 



4mJ 
Q2 • 



(3.63) 



(3.64) 



(3.65) 



3.7 The running electromagnetic coupling 

The correction due to the running QED coupling (vacuum polarization) can be implemented 

as 



(icr. 



vp 



dy 



where 



a 



a{Q') 



da 



dy 



6. 



da 



vp 



dy 



a 



A«(g2 



The correction Aa contains two parts: 



Aa{Q^) = Aai + Aa 



udcsb ; 



(3.66) 



(3.67) 



(3.68) 



the first contribution is due to the charged leptons, the second one is due to the light quarks 



{u, d, c, s, b). All details about the calculation of Sai are described in [|12|. For the calculation 
of the Aaudcsb we use the parametrization of ref. [IT^ . 



3.8 The net 0{a^) cross-section 

Now we have collected all the ingredients to construct the net QED cross- sect ion up to order 
0{a^). It has the form 



da 



QED 



dyu 



da^ ,^ , , rfa^-- 

(1 + 5vp) + 



c^y, 



, ^ Vvr"" dy^ 



, BORN , BREM ■ 

a^-^^da ^ da,^ 



+ 



dy, 



(3.69) 



M 



where the corrections S^^ are given by ([E.186|) , ( |E.196|) and (|E.215|) and the last sum in (|3.69|) 
reads 



E 



T BREM q -C/max „\/max 



^ Iv ^ ^^^ XI ^^^ ('^^ ^- + ^'^- + ^^^^^^^) ' 



(3.70) 



^ dyf, ^ JV2 

with 5"^ given by (|A.71|) , ( |A.72|) and (|A.73|) . The limits of integration in ( p.70|) are discussed 
in the Appendix B. 
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A Finite contributions to the bremsstrahlung cross-sec- 
tion 

Here we list three finite contributions, 5*^ , to the hard part of da which enter the formu- 
lae (13151 - 1071) : 
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■ Qt 
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Ami 
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e 
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1-1 
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+ 
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The indefinite integrals, in terms of which the quantities S^. , are presented here, is the first 
series of integrals over zi or 2:25 which is given in this paper. It correspontds to the first, 
innermost and the only semianalytical integration, within our numerical approach. In the 
R- frame, it corresponds to the integration over the angle ipR of the photon, see eq. ( |C.127|) . 

The three first integrals of this series are presented in Appendix D.l of |0, eqns.(D.6)-(D.8). 
Here we recall the definition of integration, and present two additional integrals. 
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The functions A, 5, C are defined in ( p.34|) -( |3.36| ). 
In the table, we introduced the abbreviation 



;max 
-1(2) 

;min 
-1(2) 



(A. 74) 



The limits of the integration z^^) ^ may be arbitrary. The particular realization of these 
limits, depending on realistic experimental cuts, which we implement within our numerical 
approach, has been derived in Appendix p.3| . 



B Numerical approach. Impleraentation of experimen- 
tal cuts 



Within our numerical approach, the cuts are applied via integration limits in ( p.70| ). There 
are cuts in variables V2, Vi and zi. 



B.l V2-dependent cut 

• Angular cut of muon 

Here we consider the cut on the angle of the scattered muon. This angle may be easily 



calculated from the Born and bremsstrahlung kinematics. In both cases 

- 2kik2 = 2 (klkl - IhlM cos9^^ 
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(B.75) 



For the case of Born kinematics one can easily derive 

.BORN S\l-y,)-2ml{Sy, + 2ml) 
cos^ = ^T- /Y- ■ (^-^^^ 

\/As\/Ai 

From eq. (|B.76|) , keeping only terms 0{m'j,), we get 






sin^ ^°— _ 1 _ ^og. ^--'- _ _L^llA_^ (B.77) 



so, one may derive the 9 . Here Vr 



BORN ,, , ,^^^ 
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v— = Sy. - '-^^^ = ^^ (l-^), (B.78) 

is the upper limit of V2 neglecting m-e. The exact expression was given in ( p.27[ ). 



Using (p.8|) , one may receive from (p.75| ) the following expression for cos 6 



,BREM 5^(1 - y^) - 2ml {Sy, -V, + 2ml) 
cosO,, = ^=^ -. B.79 

Keeping again only terms linear in m^l, we may derive the value of sine of the muon 
scattering angle in the bremsstrahlung process: 

Ami 
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sin^ ^r^" = ^^r^ (^2^^^ - ^2) • (B.80) 



It is obvious that: 

^BREM^^BORN_ ^^^^^ 

Therefore, we may introduce the upper limit of the difference 

^BORN _ ^BREM ^ ^- ^^^g^) 

or 

(C")' > (C")' - 2C"^" + d\ (B.83) 

By substitution of expressions (p.77| ) and (p.8CI| ), we receive 

T/ ^ T/c T/max \ p. ' 

V2 ^ Vo — Vo , „BORN.„ • 



/2 ^ V/g — V2 ,„BORN.„ 

^M ) 



(B.84) 



Kinematical limitation 

The kinematical limit on Vg™^^ was written in (p.27|) . 



So, the absolute max of V2 should be a minimum of the two possible maximal values of V2 

^2""^" = min [^"'^^ V2^] . (B.85) 
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B.2 Vi-dependent cuts 

• Energy recoil cut E 

The electron energy in the final state is limited from below by: E^i = p\> E , where 
E is some cut on recoil electron energy. In numerical calculation we use the SMC 
value E =35 GeV. So, 

2mepl > 2meE''''. (B.86) 

From the other hand: 

- 2pip2 = 2mepl = Ql + 2ml (B.87) 

This corresponds to 

Sy^ -V,> 2m, (E^^^ - me), (B.88) 

and we receive one of the possible cuts on Vi. 

Vi < Vr"" = Sy^ - 2m,{E'"' - m,). (B.89) 

• Energy balance cut E 

Photon energy p° also could be implicitely limited by an experimental condition, emerg- 
ing from checking of overall energy balance: 

E^ + m,-E'^-E'^=Po< E"''. (B.90) 

The SMC value is ii^ =40 GeV. It is another energy cut related to the invariant Vi. 

2meP° < 2meE''^ , (B.91) 

Vi = -2pip = 2mep° < 2meE^'' = V^"'' . (B.92) 

• Kinematical limitation 

Still we must take into account the kinematical limitation on V^^^. The boundaries of 
the allowed phase-space region are defined by ( |3.30D . 



These three cuts allow us to make the right choice for an absolute maximum of Vi 

Vr^ = min Ivr'", V^"^ , Vr""] ■ (B.93) 



B.3 zi-dependent cut 

• Cut on the electron angle The electron angle can be determined from 

- 2kip2 = 2klpl - 2\ki\\p2\ cos^e (B.94) 

It differs for the Born and bremsstrahlung kinematics. We consider first the radiative 
case, where it may be rewritten as: 

„BREM S [{Sy^ - Vi) + 2m2] - [S{1 - y^) + V2- zi] 2ml 

cos (7^ ~ 



^/X'^^fX 



{S + 2ml){Sy^-V,)+2mlz2 



%^{Sy, - V,y + imliSy, - V^] 
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:B.95) 



and 



sin^C 



Am 



AgAe 



{Ql[X,~Ql{S + ml + ml)-Z2{S + 2ml 



- mtZ2 



}.(B.96) 



In the URA in ml, it becomes 



sin^C 



Ami 



S^Ql 



S' - Qlml - S{Sy^ -V, + z^] 



(B.97) 



We may extract the correct boundaries for the invariant Zi, from three quantities: 

9g ) ' V^e ) ^iid the difference 9. 
From ( p.97| ) we get: 



/ „BREM\ 



2 Ami 



S 



S{1 - y^) + V2-z^ 
SVt, - Vi 



m. 



S 



(B.98) 



In the Born approximation Vi = ¥2 = zi = Q it reduces to 



'9. 



BORN \ 2 Am^l 



S 






m. 



S 



(B.99) 



Due to the fact that the differences between 

value, which is determined by experiment, we have: 



„BORN „BREM 

f^ — 9,_, 



BORN BREM 
9„ — 9^ 



<9. 



has to be smaller than some 



(B.lOO) 



One gets two cases 



BORN BREM 
9^ ^9^ 



(B.lOl) 



and 



BREM BORN 



With (^]98D, (gH) and ( [BJQiD we arrive at 



c\inax 



V2 + V1 






while from ( B.102 ) we receive 



c\inin 



1^2 + V^l 



(I-^m) SQl 



Vn 4m; 



9(9 + 29^ 



(B.102) 



(B.103) 



(B.104) 



Kinematical limitation 



Also the kinematic boundaries zi^^^ and 2:1™™ defined in (|3.32|) , should be taken into 
account. 
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The 9e cut conditions is active, if at least one of 

and 

(^i)"" < (^i)"" • (B.106) 

Then the absolute maximum of zi should be a minimun of the two possible maxima of zi 

z^" = min [^r", (-^i)"" ] , (B.107) 

and the absolute minimum of zi should be a maximum of the two possible minima of zi 

(B.108) 



ir" = max 



z'r, {ztr'' 



C R- frame kinematics 

In two applications in this paper, we will need another frame than the laboratory frame. When 
calculating the soft photon contribution and a table of two-fold integrals within the completely 
analytic approach, we will also make use of the so called R-frame, which is defined by 

P2+P = 0, (C.109) 



or 



Q = Pi + ki-k2 = 0. (C.llO) 

First, we introduce coordinates of 4- vectors pi, ki, k2,p, shown in the figure, 

pi = (0,0,|pl|,p?) 

ki = (O, \ki\ sinOi, \ki\ cos Oi^k^j 

p = p° {sm9Rsmipn,smeRCosipR,cos9R,l), 

P2 = {^-p'^ sin Or sin ipR,-p^ sin Or cos ipR^-p'^ cos 6 R, pfj . (C.lll) 

As may be seen from figure 6, we have choosen the z-axis of the R-frame along the vector 
Pi and matched the R-frame {z, j/)-plane with the plane spanned by vector pi and one of the 
vectors ki or k2. This trick is possible since we will use the R-frame only in analytical calcula- 
tions of tables of bremsstrahlung hard and soft integrals, in which we will always 'decouple' zi 
and Z2 using an invariant partial fraction decomposition while calculating hard integrals and 
Feynman parametrization for the soft. Finally, the photonic vector p is directed arbitrarily 
and its angles are unlimited 

0<^R< 27r, 

0<6r<7i. (C.112) 

Actually we are using two R- frames, different for zi and Z2, containing integrands with 
different pairs of 6r, ifR, both varying within the full solid angle An. 

The vector coordinates ( |C.111| ) contain, apart from angular variables, the energies Pi, k^2 
and vector moduli |pi|, |A;i2|. They can be derived using an invariant language. First we 
consider 

r = -{p2+pf = {pl+py. (C.113) 
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Figure 6: The R-frame. 



Therefore, 



pI+p' 



'T. (C.114) 

To derive an energetic coordinate of a 4-monientuni, we consider the 4-scalar product of this 
4-niomentuni with the four vector — 2(p2 + p) and write it twice: through invariants, and in 
the R-frame, using ( |U.114D . We give the simplest example: 



2{p2 + p).p = V2 — through invariants, 

= 2(p^+p°)p° = 2yrp° - in the R-frame. 



Let us introduce the following A-functions: 
Afci = \[--{p2+p-kif,r,-kl 
Afca = A -{p2+p-k2Y,T,-kl 
Api = A -{p2+p-pif,T,-pl 



[3(1 -y^)+V2f -Amir, 

\2 



Amir, 



(S - V2) 

(Sy^ + 2mlf - Amir 



A,. 



In this way, the following table had been derived: 

Afci 



l^il 
I Pi I 

\P2\ 
\P\ 



2^' 

/a7 



2v/r' 

2v^' 
IP2I, 



k\ 



"-2 



P? 



P'2 



pu 



s. 



Sil-y^) + V2 

2v^ 2v^ 



fci 



S, 



k2 



S-V2 _ 
2y/T 2v^ 

Sy^ + 2ml S:, 

2v^ 
V2 + 2ml 
2v^ ' 

V2 
2^' 



pi 



2^7' 



(C.115) 



(C.116) 



(C.117) 
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Now we write invariant S in terms of R-frame variables 

S = -2piki = 2(p°A;° - |pl||A^i| cos^i). 
By using ( |C.117|) , we get 

S 



SpiSki — y A^Afci cos^i 

2^ ' 



from which we derive an expression for cos 6i in the R-frame in terms of invariants: 

r, SpiSki — 2St 
cosyi = — . — — . 

Then we write the invariant Si 

Si = -2pik2 = 2{p\kl - \pi\\k2\ cos^2). 
By using ( |C.117|) , we get 



(C.118) 



(C.119) 



(C.120) 



^1 



SpiSk2 — Y A^Afcj cos 6^2 
2^ '' 



(C.121] 



(C.122) 



from where we derive an expression for cos 62 in R-frame: 

n SpiSk2 — 2SiT 

cos 6^2 = ^= — . (C123) 

Y Amy Afe2 

Now we can write down invariants Zi and Z2 in their two own R-frames correspondingly 

Zi = —2kip = 2p^(k^ — \ki\cos9icos9ji — \ki\sm9ism6jiCosipjij^ (C.124) 
Z2 = — 2/C2P = 2p° ( A;2 — |A;2| cos^2Cos^/i; — 1^2! sin^2sin6'/j;cosv9ijj . (C.125) 

The invariant Vi 

Vi = -2pip = 2/ (p? - \pi\ cosOr) (C.126) 

looks formally the same in both R-frames. 

Finally, relations (|C.12(: ), ( C.125 ) allow to perform / dVidz in term of / dcos9Rd(fji 



V71 



dV^ 



dzi 



ymin J^™'" TX \fRy 



V2 f + ^ f^-^ 

— d COS 6r difR. 

T J-1 Jo 



(C.127) 



This property was widely used when we discussed the treatment of the infrared divergent part. 

D Analytic approach. Tables of integrals 

D.l First and second analytic integrations. The R-integrals 



A 



L J_R 



V-," 



dVi 



dzi 
nJR. 



■.A 



(D.128) 
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In ( p.l29|) we defined the objects L^i, 1^2, Lpi, Lq^, L^, 



4m2r, 
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Ls, L^^, which will be used as 



integrands in the foregoing table. On top of them, we will also use two differences 



Dm 
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Lfc2 — L 
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D.2 The third analytic integration. The V2-integrals 
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2ymax(^^^^ + 2m2) + (V7^^)2 



T/max 

in/^°-^y^-+^\iI(t. 

(tn-l)(tii-to) ^" 
-Li2 ( —^ ] + Li2 



-Li2 



-1-tii 

ta — tl2 



t 



12 



tll)- 
ta — tl2 



In — r ln(ta - ti2, 



Lio 



-l-ti2 
ta — til 



(tl2 

- Li2 



tn-t 



11 



Li2 



■l)(tl2-to) 
ta — tii~ 

1-tll 
ta — tl2 

to ~ ti2 
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E Soft contribution to da^^ 

In ( p.52| ) the expression for da^^'^"^^ reads: 

da'""'-'' = ^^B{QlF'^^+Q,Q^F'^ + QlFlf)eie-p')dTs. (E.133) 

As ^ ' 

We substitute now the phase space dFa, (|3.40|) , which in the soft photon hmit can be 
factorized 

AT- (O \^ ^^^2 d^P2 d^p d^p f-p,nA\ 



Using definition ( p^.l4D and relation ( |3.48| ), we straightforwardly derive from ( |h].133| ) 



71 

where ( see also 0) : 

TC J 2p^ 
Then, the function F^^ takes the unique form 

^IR ^ n2 (J^ ^y , 20 [— —] [^ ^] 

^\2kip 2k2pJ '" \2kip 2k2pJ\2pip 2p2p) 

"V2pip 2p2p) ^ ' 

= QlFll + Q,Q,F'^ + QlFlf. (E.138) 



The symbol =^ in ( |E.135| ) means that instead of z-integrated IR-factors ( |3.49| )-( |33l| ), we 



use them in a completely differential form ( [E.1371 ); they depend on two photonic angles in the 
R-frame, -du, (pR, and on the photonic energy, p^. 

Since the energy of emitted soft photons is limited within a narrow interval 

0</<e (E.139) 

it is always possible to choose e small enough to ensure the phase space of soft photons be a 
sphere non-limited by experimental cuts. In other words, the phase space of soft photons is 
isotropic which allows to choose the 2;-axis along a convenient direction, differently for every 
term in ( [E.137|) . while performing angular integrations in ( [E.136|) . Therefore, the invariant 



variables Vi,V2, zi and 22 might be expressed in terms of only the polar angle cos'&r = ^ and 
made independent of the azimuthal angle ipR 

(E.140) 
(E.141) 
(E.142) 
(E.143) 
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V2 = 


-2p2p = 2p°/ = 


2mep°, 


Vi = 


-2pip = 2p\p',- 


- \pl 0, 


Zl = 


-2kip = 2p°(A;° - 


- fcilO, 


Z2 = 


-2k2P = 2p\kl - 


- hlO- 



In the R-frame 

P2 +P = 

and in soft limit p ^ R-frame degenerates to the rest frame of p2 

pI = me. 



(E.144) 



(E.145) 



This is why V2 in (|E.140|) is angular independent. 

We will use the dimensional regularization for the infrared divergences and rewrite the 
photonic phase space as follows (we took also 2p° out of every invariant variable in ( |E.14CI| )- 

(iHD) 



n 



Sp 



47r 



(20F)"r(n/2 - 1) 7o /i 



,0Nn-5^„0 i,2iT 



e (p0)n-5^pl 



,n— 4 



f^ipnTism^nT-'d^R. (E.146) 
Jo Jo 



Here 7 is the Euler constant and /z is an arbitrary parameter with a dimension of mass. 
Now integration in (pi).136| ) over p^ is performed straightforwardly: 



+1 



s:IR,soft 



d^ 



Pi^ + lni + iln(l-a 



T' 



IR 



In ( Fl4^ ) 



and the typical pole term 



0n2 77IR 



5IR 



T'"^ = {2pyF 



11,1 

h -7 + In — ^ 

n-4 2' 20F 



(E.147) 



(E.148) 



(E.149) 



represents the infrared divergences at n = 4. 

From ( |E.14CI| )-( |E.143| ) entering ( |E.138| ) and from ( |E.147| ) we see that only those integrals 
over ^ may occur, which are presented in Appendix D.2 of 0). Writing JF^^ in a form similar 
to (|E.138|) , we will calculate the three contributions (5|,^'''°^\ Sj^'^°^^ and 5^^'^°^* separately. Before 
this, however, we present a collection of formulae in the R-frame in the soft photon limit. 



E.l R-frame kinematics for the calculation of the soft photon con- 
tribution 



The R-frame is defined by 



P2 +p = 



(E.150) 



or 



Q = Pi + ki - k2 = 0. 



(E.151) 
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Since the R-fraine is isotropic, there is no need to fix its 2;-axis along a given direction, say 
along pi0 



Pi 



0,OApIIPi 



(E.152) 



It might be equally chosen along ki, 

h = (0,0,\ki\,k',), 

or along k2, then 

k2 = (o,o,|fc;u°), 

— # — * — * 

or along any linear combination of any vectors, say ka = a^i + (1 — a)k2, then 

fc„ = (0, 0, \aki + (1 - a)k2\,ak^^ + (1 - a)^;") • (E.155) 

So, we have indeed many R-frames, which differ one from another by a spatial rotation and 
when we write an arbitrarily oriented photonic 4-momentum as 



(E.153) 



(E.154) 



P 



p°(sin ^i? sin ipR, sin {^r cos ipn, cos {}r, 1), 



(E.156) 



one should understand that in every R-frame one has its own angles "i^ij, ifR which vary within 
the same limits - covering the full solid angle. In this way, we arrive at equations ( |E.14CI| )- 
( Pj^.143 ) for invariants V2, Vi, Zi, Z2 with formally one parameter ^. 

In the expression ( [E.147D , which has to be integrated over ^ with ( |E.148| ) and ( pi], 137] ), 
enter the energies p?, k^2 and moduli |pi|, 1^12! ( see ( [E.140|) -( [E.143|) ). 

All the energies and momenta moduli depend only on three invariants: S,y and V2. We 
make no distinction between y and y^ in the soft photon kinematics. In the soft photon 
problem, we neglect the small invariant V2 as compared to the others. In this limit, the 
table (ICllTp reduces to 



l^il 



2me' 



2mp 



ki 



rvn 



Si 
2m f.^ 



S 
2mp 



A? 



\fl\ - 


_ V "e 

2me ' 


P? 


IP2I = 


V2 
2^7' 


Pi 


\p\ - 


= P2 , 


po 



Sy + 2m; 



2r7i„ 



rrie. 



2mp 



^In this subsection all 4-momentum coordinates are understood in the R-frame. 



(E.157) 
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E.2 Muonic current 

In muonic current we are dealing with 



T-IR 



h 



k. 



2kik2 



2kipJ \2k2pJ {2kip){2k2p)' 
see ( pi].137D . By applying the Feynnian parameterization for the last term, we have 

2kik2 



(E.158) 



7IR 
MM 



mz 



mi 



da- 



(E.159) 



{2kipy {2k2pY Jo {2k^pf 
where a new 4- vector was introduced as discussed in (|E.155|) : 

ka = aki + (1 - a)k2. (E.160) 

Therefore, using all soft machinery described above we can write 5^^*^'^°^* as follows 



rIR,soft 

MM 



Qj|(pi^ + ln£)i/rfa/rfe^J + i/rfa/rfein(l-a-^?}, 



see ( F147D , with 



-F, 



IR 



m. 



m. 



Sy + 2m 



kfii-m' kri^-m' 



UO' 



[i-Pciy 



On passing, we used 



Z2 



2k^p = 2/(A;°-|4|0 = 2pXil-PaO, 
-2k,p = 2p%k',-\h\0 = 2p0fc?(l-/?iO, 
-2k2P = 2p%k'2-\k2\0 = 2/A;°(l-/?20, 



with three velocities 





fell 


VA/ 


/3i = 


k\ 


s, ' 


P2 = 


H 


V^s 


k\ 


s 


p, - 


\ko\ 





1,0 



(E.161) 



(E.162) 



(E.163) 
(E.164) 
(E.165) 



(E.166) 

(E.167) 
(E.168) 



We also have from ( [E.16C ) the following relations: 



fc^ 



Sa + 5*1(1 — a) 

2mf, 
m? + a{l — a)Sy. 



(E.169) 
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Therefore 



/So 



[Sa + 5*1(1 - a)f - Ami ^l + "(1 - (^)Sy 



Sa + 5'i(l — a) 
Using the table of integrals from Appendix D.2 of [^ we obtain 



(E.170) 



rIR,soft 
MM 



Qt da{ P'^ + \n 



2e 



m. 



m. 



1 m. 



In 






fcf (1 - P!) kf (1 - /3: 



+ 



Sy + 2m 



k^' 



:i-/5; 



1 m. 



In 



1+/32 



2/3i fco^ (1 - PI) 1 - /5i 2/^2 A;0' (1 - Pi) 1 - /32 
1 Sy + 2ml 1 



In 



l + /^a 

2/3„ fcf (1 - /32) "' 1 - /5„ 
We have from (|E:T66|) - (^168D 

{k',)\i-P!) = 



(K 



0\2/ 



A 



-kf 

— Kg 






And the expression ([E.161|) becomes 



rlRjSoft 
Vm 



H — TT In ;^ + ^r-^ In ■ 



c/a 



m^ + «(! — q;)S'?/ 



2A l-A 2/^2 1-/32 
1 



+ 



S'y + 2m 



(ia 



In- 



/?a 



2 7 /?,[m2+a(l-a)^y] 1 + /?„ 
Now we use the URA in the electron mass 



1 , 1+A 
In ^ 

2/?i 1 - Pi 

1 , I+P2 
In ^ 

2/32 1 - P2 
The first integral we calculated precisely in | 



In 



^1 



m^m^ 



In 



S 



m^m^ 



(Sy + 2m; 



(ia; 



[m; + a(l - a)^?/] 



l + /3^ 

/3 



with 



P 



\ 



In 



Am,f, 
1 H ^ 

/3 + 1 



/? 



(E.171) 



(E.172) 
(E.173) 
(E.174) 



(E.175) 

(E.176) 
(E.177) 



(E.178) 



(E.179) 
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With the second integral the situation is more comphcated. Defining it as in 



bs, = 



1 )■ 

- {Sy + 2ml) J 



da 



In 



I -Pa 



f3a[ml + a{l-a)Sy] 1 + (3^ 



(E.180) 



we cannot take over the analogous result from since there that integral was calculated in 
the URA in the leptonic mass m. Here we need result exact in m^^ (see also 0). Making use 



of URA in nie, with the aid of (|E.17CI| ) we get 

1 
S^ ^ (^Sy + 2m 



da 



^^7 ml + ail-a)Sy 



In 



m; 



m^ + a(l — a)Sy 



[Sa + Si{l-a)Y 
The expression for 5*$ simplifies drastically and can be calculated straightforwardly 



(E.181) 



5, 



Sy + 2ml( 



m?{Sy + 4m^ 



In^ 



2 
+Li2 



52(1 - y){l - yai){l - ya2 

012 



In 



^2 



-ttlj 



(-«l)(l-y) 

;-ai) 



+ ln(l-y)lnii^^lM^-Li2(ll 



;i -ya2) 



y)oi2 



with 



+ Li2 



"1,2 



'l-y){-ax) 

Ci2 



lT/5 



\An 



-OiV 

02 / 



Collecting all terms together, we finally have 



rlRjSoft 
MM 



Qj 2 Pi^ + ln 



L^-1 



In 



^'(l-y) 



"^e"^^ 



5, 



# 



/i/ V 2/3 
From the virtual photon correction to the muon vertex we have the contribution 

m, 



(E.182) 



(E.183) 



(E.184) 



rvert 
MM 



IR 



2(3 



Lp In 75:5 7 + L12 I :; ) - L12 



q^-2P^^ + ln-^ 



P^ 



1-/5. 



1 + /3, 



(E.185) 



The infrared divergence and the scale parameter /i cancel exactly in the sum of these two 
contributions. The sum reads 



'/i/i 



IRjSoft I xvert 

'2 



^iR,sott ^ ^ 



'M 



1+/?^ 



2/3 



' T^min \ o 

21n^-ln(l-|/)j+-/?L^-2 
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L^ln- 
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1 



i/tti)(i -ya2) 



\n\l - y) 



+ ln(l-,)lnil^#)il:i^ + Li2 



'1 - 2/02) 
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(1 -2/)tt2 



+Li2 



■(l-Z/)(-«i)' 



2Li2 



02 



02 



(E.186) 
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E.3 Electronic current 



The correction 5^^'^°^* is very easy to calculate. 






+ 



P2 



2piP2 



2p2pJ {2pip){2p2p)' 



(E.187) 



Since 2p2P is independent of ^ = cos i)^, the a-parameterization is not needed here and S^^'^°^^ 
reads 



rlRjSoft 
ee 



P'"" + In ^) ^ Jd^:Flf + \ Jdi ln(l - e):F\ 



Introducing 



we receive 



\pi\ p'y' + ^mlSy 
^ -0 Sy + 2mi ' 



2^^^IR 

ee 



IR 

ee 



m; 



Pi 



1 + 



(E.188) 



(E.189) 



2mePi 



pf (1-/30' rn^P'i (l-/?0 



■1 + 



r/i' 



i-/3e pf(i-/?0'' 



Using the table of integrals from Appendix D.2 of M we derive 



s:IR,soft 



Q^ P- + In^ 



-l + iln.^ + ^ ^' ' 



/3 l-(3 pf(l-/?2) 



(E.190) 
(E.191) 



+^ + ^ 



Li2 



' 2/3 ^ 



Li2 



^ 2/3 - 



+ 



m? 1 



Pi 



o2 



1 .i„.i + '' 



2/3(l-/J2) 1-/3 



In the URA in m^ we obtain 



and 



Li-; 



' 2/3 ^ 



1 1 + /3 , 5y 
- In ^ 2 In ^ 

/3 1 - /? m2 

■ 2/3 ■ 



- Li2 



/3+1 



^-2Li2(l)-21n2^. 
mi 



(E.192) 



(E.193) 



^iR,soft ^ 2Q^<|(pi^ + ln-) (-l + ln^)+l + ln^-ln2^-Li2(i: 



So, we can write 

2e\ f ^ , Sj/ 

(E.194) 
The corresponding virtual photon correction to the electron vertex has the following form: 

(E.195) 
The complete answer is: 



Q 



IR,soft I rvert 

Sy 



^lK,sott _^ ^ 
"ee ' "ee 



In^lf-l) (ln^ + 21n^ 
mi J \ mi mi ^ 



3 , Sy I ^ n Sy 
1 + - In ^ - - In^ ^ 



2 m^ 2 



mt 



. (E.196) 
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E.4 lie interference 

In the /ie interference the expression for F^ reads 

2kipi '2kiP2 2/C2P1 



AR 



+ 



2/C2P2 



^' {2kip){2pip) {2kip){2p2p) {2k2p){2pip) {2k2p){2p2p)' 
We introduce the a parameterization with the aid of two new 4- vectors 

kia = kia + pi{l- a), 
k2a = k2a+p2{l- a), 

resulting in 



7IR 

fie 



Si 



s 



2A;V(l-/3iO 2k^2P'{l-P20 



+ da 



Si 



S 



2{k2apy 2{ki^py 



(E.197) 



(E.198) 
(E.199) 



(E.200) 



and S]j^'^°^^ becomes 



xIR,soft 
fj,e 



QuQ 



fi^e 



1 +1 



1 +1 



P'"" + In ^j^ da J d^J^'J^ +y da I d^Hl ~ e):F'^ 
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fie 
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IR 
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-1 



^1 



i-M i-M {ki:)\i - p2aiY {k\:)\i-fiuy 



(E.201) 



(E.202) 



-k{^ = mla + ml{l - a) - 2kipia{l - a) 



'^Ic 



mf^a + mlil - a) + Sa{l - a) 
Sia + {Sy + 2ml){l-a) 



2raf. 



(E.203) 



—k2a = rnj^a + mKl — a) — 2k2Pia{l — a) 
= rri^^a^ + rnlil — aY + Sia{l — a), 
Sa + {Sy + 2ml){l-a) 



k^ 



2m^ 



(E.204) 



P: 



la 



I "'la I 



(E.205) 



Sia + {Sy + 2ml){l-a) 



Ami 



mla^ + mlil - af + Sa{l - a) 



Sia + {Sy + 2mi){l 



a) 



A 



kic 



'2a 



1,0 
'^lo 



(E.206) 



Sa+{Sy + 2ml){l-a] 



Ami 



m^a^ + ml{l - af + Sia{l - a) 



Sa + {Sy + 2ml){l-a) 
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Using the table of integrals from Appendix D.2 of [^ we derive 



rlR.soft 
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(E.207) 



S 



Si 



1 (ia 

^F" 



la 



5 



(ia 



In 



1 (ia 

— /i;2a 
l-Pla 



2 -^ ^^la/^la 1 + /?la 

5 /■ da , 1 - /52a 
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(E.208) 
(E.209) 
(E.210) 
(E.211) 



Where we introduced the generic function S'$(/, /) 
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2 J mla"^ + ml{l - af + Ja(l - a) 



In- 



m: 



m^a^ + mlil — a)^ + Ia{\ — a) 



/a+(5y + 2m2)(l-a) 



Mn^lu '" 



In 



2"^; 



om- 



2^2 A — I 



2 S'2?/2 mlm^^ 4 

+ ln — ln-^--ln2— . 
5?/ / 2 ^1/ 



jl„^-^-LMl 






The expression (|E.207|) reduces in the URA in rrie to 



(E.212) 



rIR,soft 
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In the next Appendix on the box contribution we derive 



rbox 



QeQ 



C^fl 



2 /i2 



Pi^-^ln^ 41n(l-y) 



+ B 



fin 



and the following short final expression is obtained 



(E.213) 



(E.214) 



K 



Ke = 5" + C^ = QeQ, 41n(l-y)ln^ + P^ 



Sy 



fin 



(E.215) 



The finite contribution from box diagram, Br;^, is presented in Appendix F. 
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F Two photon exchange contribution 

The two-photon exchange contribution is described by the two box diagrams: direct box and 
crossed box. For the sake of symmetry, it is convenient to duphcate the number of diagrams 
and to deal with four diagrams shown in the figure below: 



h 



a 



ki+P 13 ^2 



ki (3 k2-p 



a 



k. 
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q + p 



q + p 



P 
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h+P (3 k2 
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ki (3 k2-p 



a ^2 




Pi-P 



This is why the two-photon exchange matrix element is given by an expression, which 
contains an extra factor 1/2: 



M^ 



d*p 



2 J p^ij) + qY 
xu{k2 



'yf3{2kia + Pla) {'2k2a - laPh/S 



p^ + 2pki 



Xu{p2 



(2p2a + laPhp l^Pla 



p"^ — 2pk2 

PI a) 



p^ + 2pp2 p — 2ppi 

For later use we introduce the short hand notation for propagators 



u{ki) 
u{pi). 



IK^i) 


= p^ + 2pki, 


im 


= p^ - 2pk2, 


n(pi) 


= p^-2ppi, 


life) 


= p^ + 2pp2. 



(F.216) 



(F.217) 

Now we separate the infrared divergent (IRD) part of the two-photon exchange contribu- 
tion. There are two IRD-poles: 1) one is located at p — > and 2) another one - at {p + q) — * 0. 
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At p -^ 0, one has 



1 1 



2n2 



2p2g 



4fciP2 



Akipi 



4A;2P2 



Ak2Pi 
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7/3® 7/3, (F.218) 



While at p — * g, one has 
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7/3 ® ll3- 



(F.219) 



For the latter term, we perform the substitution p = p —q and observe that at this substitution 
the propagators transform as follows 



(^l) - 


- n(-^2), 


ih) - 


- Ei-h), 


iPi) - 


- n(-P2), 


iP2) - 


- n(-pi)- 



(F.220) 

With one more substitution p -^ — p, and using equalities /c2-P2 = ^i-Pi, ^2 -Pi = ki.p2, we see 
that the second pole gives exactly the same contribution as the first one. Therefore, the full 
IRD-part of the two-photon exchange reads: 
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7/3 7/3. (F.221) 



Now we add and subtract the IRD-part of the two-photon exchange and when adding it 
we use ( [l''.221| ), while when subtracting it we use the sum of ( |1:''.218| ) and ( [1:''.219[ ). Using the 
identity 



1 
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p^{p + qy 



p2q2 



1 _ -2p.{p + q) 

q^ip + qY q^p'^{p + qY 



(F.222) 



we arrive at the final expression for the two-photon exchange contribution before integration 
over df'p. 
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Apjp + q) r -S 

p'^ip + qYq'^ 

1 r kia 



7/3® 7/3 
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n(^i)n(pi) n(^i)n(P2 



+ 



+ 



'^2a 



p'ip+qriuih) 11(^2) 



p2(p + g)2 
1 



ll3Pla laPlp 



IpPla laPlp 



7/3® 
®7/3 



2p'^{p + g)" 



IK^i) 11(^2 



7aP7/3 

[\{{P2) 

_P2o_ 

\{{P2) 

laPlp 



+ 



7/3® 7/3 

n(pi). 

Pla 

n(pi). 



\{{P2) \{{Pl 



(F.223) 
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The first raw of tliis formula represents the IRD-contribution, the second raw stands for 
the IR-free scalar, the next two represent the IR-free vector, and finally the last one is an 
IR-free tensor. 

For the IR-divergent part we introduce the well-known presentation 

167r^ f d'^p 



J3 {ki,Pi) = — / 



(27r)Vn(A;i)n(F 



Jo {-K^} 2 Jo {-K^) q^ 

where 

-Kl = my + ml{l-yf-Sy{l-y). (F.225) 

We will use the short hand notations 

dy 



pIR ^ pIR(^) ^ _p^^(^) f' 

^ u 



{-Kiy 



Ks = KiS) = /-^ln^4^, (F.226) 



Then all needed IRD-integrals are 



,IR , Ks 



■h\kuPi) = J3 (^2,P2) = -Pr + ^, (F-227) 

jT{kuP2) = jTik2,Pi) = P^ + ^. (F.228) 

with 

Ks, = K{-Si). (F.229) 

For the IR-finite scalars 

we derived 

JA{ki,Pi) = .Mk2,P2) = Ks, (F.231) 

J4{ki,p2) = J4{k2,Pi) = Ks„ (F.232) 

F.l Vectorial integrals 

We begin with the reduction of vector integrals. First we define the coefficients of the vector 
reduction 



V. 



ij _ 1^^ f d^PPtJi 



% J {2nYp'{p + qyU{h)U{Pj) 



{v,rq, + {Vkritd, + (v.ripj),, (F.233) 
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The system of linear equations for these coefficients in terms of scalar integrals, derived by 
contraction of ( |1''.233| ) with g^, (fcj)^, {pj)f_i, reads 

2gV;^' + 2Sk,q ■ k,V^' + 2S^^q ■ p,V;^ - Mh,p,) + 2f^{h,p,) = 

2Sk,q ■ ky^ + 2Su;'i^vi^ + 2Sk,s,^h ■ PjV;^ - Up,) + J7\h.P,) = o, 

2S,^q ■ p^^ + 2Sk^S,^h ■ PjVl^ + 2S,^^p]V;^ - Uh) + J^ih^P,) = 0. (F.234) 
Here 

Ski = ~^} Sk2 = +1, 

Sp, = +1, Sp, = -1. (F.235) 

Two additional types of scalar integrals 

167r^ f d^p 

'^ ('•' = — I (2x)v(P + «)^n(fe)' '^-'''^ 

and 

IGvr /■ dp 

'' ^p^^ = -r J {2n)v{p+,m{p,y ^^-'"'^ 

are introduced. They safisfy the equalities 

^3 (ki) = h {k2) (F.238) 

and 

J3iPi) = J3iP2), (F.239) 

and may be described by one generic formula: 

J ( \ = [^ '^ i^ (g ~ ii)'^yi^ -y)- qK^ - y) 

Jo -q^y - gf (1 -y) + {q- qiYy{l - y) q^y 

(F.240) 

with qi = {ki, -k2, P2, -pi)- 

The solution of the system (|1^'.234D is presented in subsection F.4 of this Appendix. 



F.2 Tensor ial integrals 

Define the coefficients of tensor reduction 

167r^ /■ d^pp^py 



rplj 



^ I i2n)^p'ip + qyUiki)UiPj] 



= TqS^^ + Tglq^q^ + Tli{ki)^{kj)^ + T'^^p{pi)^,iPj)u 
+Tqk [%{ki)v + q^.{k^)^] + T^^ [q^{j)j)u + q^.{Pj)^ 
+Tl[{k,),{p,\ + {kUPi),]- (F-241) 

The system of equations for the coefficients T*-^ was derived in the following way: 
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the first equation was received by multiplying ([F.241|) with S^,^ 

the next three equations were received by multiplying ( |1''.241| ) with q^ 
d^p [2pq + q^ + p^ — q^ — p^] p^ 



(F.242) 



{2TiYp\p + qYmki)Y{{P3 
d^p p^^ 



{2n^)p'U{h)U{Pj 



d^p 



pp 



{27r^) {p + qyU{h)U{Pj] 



d^p 



p^ 



{2'K^)p\p + qYl{{h)Y{{p,] 



= 2qX' + '^q'q.Tg + q'Tiiik^, + q'r^^{p,), 

+T% kg2 ^ 2q\h),\ + Tg \q,q^ + 2q\p 



^j/M 



+g2Tg [{h), + {p,),] , (F.243) 

• another six equations were received by multiplying ([F.241| ) with (/cj)^ and {pj)u- 
The latter nine equations are 

o'Tiil + 2q'T% + q'Tg + gV,(fc„p,) = 

2kJTg + q'T^ + 2Sk^S,^h-p,Tg + V,,{h,p^) + gV,(fc„p,.) + Jr(^.,P.) = 

2S,^S,^h ■ PjTii + q^T% + 2p]Tl - V^uih^p,) + V,,{h,Pi) = 

q'T^ + 2g2T^^ + g^Tg + qX{h,p,) = 

25fc,^5p^.fc, ■ p,T;j, + q'Tg + 2klTl - Vspih,p,) + V,pih,p,) = 

2p^ ■ pfl'i^ + g^T^i + 2^fc,^p^i;, ■ p,T| + V^ip(A;„p,) + Js'^lA;.,^,) + q^V,{h,p^) = 

2g2T^^^ + (fT% + g^T^^ _ V^k{h,p,) - Vi,{h,p^) + 2gV,(A;„p,) = 

q ^qq + '^^i^qk + '^^ki^pjh ■ PjT^p 

-V3qiki,pj) - Vik{ki,pj) - Vip{ki,pj) - J^^iki.pj) = 

qXi + ^S,^Sp^h-p,T;i + 2p,Xi, 

-V2q{ki,pj) - Vik{ki,pj) - Vip{ki,pj) - f^{ki,pj) = (F.244) 



Together with the tenth equation ([F.242| ), they can be solved and the answer may be expressed 
in terms V^^ and additional vectors arising from various pinches of the box diagrams. 



F.3 The box contribution to the cross-section 

Substituting all these solutions to the box amplitude and computing its interference with 
Born amplitude, we derive the contributions to the unpolarized and polarized cross-section 
parths from the two-photon exchange diagrams to the differential cross-section of elastic fie 
scattering 



da. 



BOX 



dy 



-iP^^in) ln(l - y) 



da 



BORN 



dy 



~r l^unp ~r -i^e-^fj.'^pol 



(F.245) 
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where unpolarized and polarized non-factorized contributions Bunp and Bpoi are: 

2 „ / 1 \ 



i3. 



with 



unp 



y y 



y 



2i?^ 1 



y^ii 



\y J \y Vi^J 



1//.1 



+4i?. 



+2 



1 I 2(2-y) 



M/j 






S{Ks - yJp) 



\y y + ^Ra, 



SJk, 



E^ 



'pol 



+4 
+2 



-1 + --2RA{SKs + S,Ks,] 

y J 

1-R 



-1 + 2R, 



y 



-{1 + 2R, 



-1 + - 
y 



2Rl 



V- 



Is + 2 



— Rn 



+4fl„fl + i^^ 2<l + 2«"' 



y^^l 



t/ + 4i?. 



2i?^(l - R^) 
- i?^ In ^ 



1-2R, 

y 



^Si 



+4 



R^yjl + 2i?, 



-SJ, 



K, 



(F.246) 



^ S ' 

1 



1 + i?. 



There are two useful combinations: 



SiKs, + SKs 
S{Ks - yJp) 



-In 






2Li2 1 



m. 



S 



Q 



Q 



m. 



TT 



In" ^ + - In" ^ + 2Li2 1 + ^ + ^• 



S 



m. 



S 



n 



(F.247) 

(F.248) 
(F.249) 



Here 



and 



Jk 
Jp 



J{m 
J{m\ 






m 



In— I In 1/^2 -2112(1-1/^2) --In y^.,+Ti 



(F.250) 
(F.251) 





Q' + 2m2 + , 


J^l 


y^l2 = 


2ml 
1 





(F.252) 
(F.253) 
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The two-photon exchange contribution simphfies drastically in the DOURA. We list the 
answer because of its elegance. 



, BOX 2 J- 



dy S 

+ (1 - PeP,) 



4 + (l + PeP,)(l-- 

y V y, 



81n(l-y) 
(l-'^j{\n\l-y)-2\nil-y)\ny + 



vr^ 



2 
+ -ln(l -y) -2\ny 

y 

+2 [1- -J \n^{l-y) -4^^1n(l-y)l. (F.254) 

G Description of www-Figures 

In March 1996, when the FORTRAN code fiela was completed, a lot of figures for the radia- 
tive corrections to the polarized cross-sections and the polarization asymmetry were produced 
and put to a home-page of the theory group of DESY-Zeuthen, pittp : //www . if h . de/[ ^bardin. 

In this Appendix, we give a short guide to these figures. 

All the cross-sections, radiative corrections and asymmetries are shown as function of y^. 
The units are: cross-sections in microbarns, radiative corrections in percent, asymmetries in 
absolute, dimensionless units varying from to 1. 
Some legends: 

- parallel, antiparallel means mutual longitudinal polarization orientations; 

- Pe and Pfj, are modules of polarization degrees, always written in figures; 

- All corrections mean that all 12 contributions ( |1.4|) are taken into account; 

- AN means Analytic calculations; 

- when AN is not written, the numerical calculations are meant. 
The 6y^ is defined by 



^.= 


, QED 

da 
dyi. 


, BORN 

da 



1, (%). (G.255) 

dy^^ 
SIURA-series of 43 figures contain following plots: 

• 1-3 Born cross-section and Born asymmetry; 

• 4-5 Results of completely integrated analytic calculations without cuts; 

• 6-7 Results of completely numerical calculations without cuts; 



53 



• 8-17 Illustrate effects of some cuts separately, muon angular cut being treated both 
analytically and numerically; 

• 18-19 are our main result: All corrections, All cuts (four indeed); 

• 20-27 illustrate four cases: fi + e corrections, fi corrections, e corrections and /xe inter- 
ference corrections, correspondingly, without cuts; 

• 28-35 the same, but with all cuts; 

• 36-43 is a series of figures with 'realistic' Pg and P^. It is quite senceless for the asymme- 
try, where Pg and P^ cancel. We only 'gain' an instability due to cancellation of nearly 
equal numbers. 
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